Noise-correlation-time-mediated localization in random nonlinear dynamical systems 
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We investigate the behavior of the residence times density function for different nonlinear dy- 
namical systems with limit cycle behavior and perturbed parametrically with a colored noise. We 
present evidence that underlying the stochastic resonancelike behavior with the noise correlation 
time, there is an effect of optimal localization of the system trajectories in the phase space. This 
phenomenon is observed in systems with different nonlinearities, suggesting a degree of universality. 
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By stochastic resonance (SR) it is normally understood 
the phenomenon by which an additive noise (usually con- 
sidered uncorrelated) can enhance the coherent response 
of a periodically driven system. First proposed in climate 
model studies [Q , SR was first experimentally verified by 
Fauve and Heslot @, and since then this behavior has 
been predicted and observed in many different theoret- 
ical and experimental systems (see || for an extensive 
review and a complete list of references). In particular, 
the presence of SR has been discussed in a great number 
of models including spatio-temporal systems Q , and has 
helped to understand how biological organisms may use 
noise to enhance the transmission of weak signals through 
nervous systems [j5|,^| . Quite recently it has been numer- 
ically shown that SR can also occur in the absence of 
an external periodic force as a consequence of the intrin- 
sic dynamics of the nonlinear system Q , a behavior that 
has been denominated autonomous stochastic resonance. 
Most of the work on SR has traditionally focused on sys- 
tems with additive noise, and with some exceptions (see, 
for instance, Q) little attention has been given to cases 
where the noise perturbs the system parametrically, in 
spite of the well known differences with the additive sit- 
uation. With respect to non-white noise, the effect of 
additive colored noise on SR has been considered in peri- 
odically driven overdamped systems |J , showing that the 
correlation time can suppress SR monotonically, a fea- 
ture demonstrated experimentally in However, only 
very recently the situation in which the system is sub- 
ject to both multiplicative and colored noise has been 
discussed in the literature. In jll| the authors analyze 
the effect of multiplicative colored noise on periodically 
driven linear systems, discussing the appearance of SR 
by changing either the intensity or the correlation time 
of the noise. For nonlinear models, in jl^] we consid- 
ered a system without periodic external force but with 
an intrinsic limit cycle behavior, which was parametri- 
cally perturbed by an Ornstein-Uhlenbeck (OU) noise, 
finding a nonmonotonic behavior of the coherence in the 
system response when measured as a function of the noise 
correlation time, while no coherence enhancement was 



obtained when changing the noise intensity. A similar 
result has also been recently obtained analytically for an 
overdamped linear system periodically driven and para- 
metrically perturbed by an OU process ]l3]j . 

In this paper, we present numerical evidence which 
suggests that underlying the SR-like behavior as a func- 
tion of the noise correlation time, there is a localization 
effect of the system trajectories in the phase space for a 
particular value of the correlation time. This is obtained 
in systems with intrinsic limit cycle, perturbed paramet- 
rically by an OU process and with different nonlinearities, 
which is also a clear indication that the phenomenon is 
not a peculiarity of an specific model. 

We study three different 2D random systems. The de- 
layed regulation model, known from population dynamics 

s 

x t +i = A t a; t (l - x t -\), (I) 

the Sel'kov model for glycolysis |l5| 

x = —x + \ t y + x 2 y, 

y = b-X t y- x 2 y, (2) 

and the Odell model also from population dynamics |llf | 

x = x[x{\ — x) — y], 

y = y(x-\ t ). (3) 












FIG. 1. Phase space partition of Eq. (|), (a) with N+l = 9 
deterministic periodic attractors and (b) with superimposed 
random states (dots), and with ay = 0.05 and r = 3. The 
thick line corresponds to the attractor T((A}). 
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FIG. 2. Residence times density function for (|2j), obtained 
with 50 realizations of 5 x 10 7 time steps of size At = 10 -2 , 
<A) = 0.123, a = 5 x 10~ 4 and r = 6. 




FIG. 3. RTDF versus r for (a) Eq. (@), (b) Eq. (|), and 
(c) Eq. (|). Data obtained with (a) (A) = 2.02, oy = 0.05 and 
100 realizations of 10 6 iterations, (b) (A) = 0.123, a = 5xl0~ 4 
and 50 realizations of 5 x 10 7 time steps of size At = 10~ 2 , 
and (c) (A) = 0.49, a — 10~ 3 and realizations as in (b). The 
inset plots indicate the height and width behavior. 



Here, t takes discrete values in (gj) or continuous values 
in (||) and (0), and in all cases we will consider the con- 
trol parameter as a random variable A t = X + Ct, i- e -, as a 
deterministic part A, plus an stochastic perturbation £ t , 
which is assumed to be an OU process, i.e., a stationary 
Gaussian Markov noise with zero mean, = 0, and ex- 
ponential correlation (CtCtf) = (D/t) exp (— \t — t'\ /r), 
where r is the correlation time and D/t = a 2 is the 
variance of the noise. We will refer to the square root 
of the variance, tr, as the intensity of the noise. The 
deterministic counterparts of (|l|), (|2|) and (|3|) undergo a 
supercritical Hopf bifurcation at A = Xh which, in the 
Sel'kov model, also depends on the parameter b. 

The numerical integration has been carried out with 
A in the limit cycle parameter domain. The iteration of 
(|l|) has been recreated using an integral algorithm |l7|] 
that guarantees the quality of the correlation function in 
the simulations of the noise at discrete times, while (||) 
and (||) have been integrated by an order 2 explicit weak 
scheme . The results presented hereafter are indepen- 
dent of the initial conditions and were obtained after the 
decay of the initial transients. 

The observed fact |l2| that for a particular correlation 
time r r the coherence of the system oscillations has a 
maximum, and that the frequency of these oscillations is 
close to the deterministic one Wd, seems to indicate that 
for the resonant correlation time, the probability that 
the system visits the attractors associated with the mean 
control parameter value A =< X t > has also a maximum. 
If this is the case, this maximum should be accompanied 
with a decrease in the probability to visit other attrac- 
tors associated with parameters far away from A, or, in 
other words, should lead to an effect of concentration 
or localization of orbits around the attractor associated 
with A as soon as t ~ r r . It is worth recalling that be- 
cause of changes in the stability properties, a somehow 
similar localization effect can also occur in parametric 
deterministic systems with time dependent parameters, 
as is the case, for instance, in the well known parametric 
resonance phenomenon. 

With the aim of studying the residence times distri- 
bution of the system on the different available attrac- 
tors in the system periodic or quasiperiodic domain, 
we consider a deterministic attractor A(A), i.e., the at- 
tractor obtained with the deterministic counterpart of 
the stochastic system, evaluated at a particular value of 
the control parameter, A. Next we divide the system 
phase space in N + 1 attractors associated with N + 1 
values of the parameter separated a distance AA. In 
this way, a mesh is composed by concentric determin- 
istic attractors centered around the stationary equilib- 
rium state (x*,y*) |A~Aff7 with A in the fixed point do- 
main. This partition looks like the one shown in Fig. 
la. With this construction, we have a series of N + 1 at- 
tractors {A(A JV /2_)-A(Ai_),A(A ),A(Ai + )...A(A iV/ 2 + )}, 
where we use the definition Afe± = A ± kAX. This series 



divides the phase space in N rings, each one denoted by 
r(7 fc ) = (A(A fe ), A(A fc+ i)), where j k = (X k +i + A*) /2 is 
the mean control parameter obtained with the control 
parameters that define the boundary of the ring. The 
stochastic system is integrated on this mesh, and its evo- 
lution describes random trajectories as the one described 
in Fig. lb for the particular case of (Q), visiting during a 
finite time each ring of the mesh. During the integration 
process we measure the residence time in the rings as fol- 
lows: let t\ and t\ be the entrance and exit times to the 
ring r(7fe), respectively. The residence time in this ring 
is i(7fc) = t\ \ — t\ , and we denote the residence time of the 
n visit event to the ring r(7fc) by fn(7fe)- Then, if during 
an integration time /, which is achieved by integrating 
R realizations of M time steps, there have been 14 visit 
events to the ring r(7fc), the mean residence time of the 
system in this ring is given by the mean of the residence 
events, that is, T(T(j k )) = Y%Li H^- Such a determi- 
nation of the residence times gives an alternative statis- 
tical measure of the resonant amplification described in 
|^2| . Therefore, given a pair (cr, r), the function defined 
by the histogram P(T) = P(T( 7fe )) = P( T( ^ fc)) ) is a 
measure of the probability density for the system state 
to be in the region defined by the ring F(7fe) . An ex- 
ample of histogram is depicted in Fig. 2 and shows that 
the system mostly visits the attractors surrounding the 
ring r(< At >). We remark that we have carefully se- 
lected the simulation parameters to ensure that the par- 
tition does not contain overlapped attractors such that 
this has a well defined meaning. An illustrative example 
of the residence times density function (RTDF) as a func- 
tion of the correlation time is depicted in Fig. 3 for the 
three models. Obviously, the localization of the system 
trajectories depends strongly on r. The RTDF height 
shows a nonmonotonous behavior reaching a maximum 
at a particular value of r ~ r* and, at the same value, the 
width W calculated at the height h/y/e shows a remark- 
able minimum, as represented in the inset curves. The 
correlation time of the parametric random perturbation 
acts as a tuner which controls (in an statistical sense) 
the behavior of the system, maximizing its localization 
on the region of the phase space surrounding < At >. 
Furthermore, the relation h/W has a maximum for a 
particular value of r, and this optimal value depends on 
A =< A t >, as can be appreciated in Fig. 4. Such a depen- 
dence enables us to relate the optimal correlation time for 
maximal localization, r*, with the temporal scales of the 
deterministic counterparts. We first study the behavior 
of the postponement of the bifurcation point because of 
the multiplicative noise in order to obtain the postponed 
bifurcation point A*y(<7, r).We next calculate the effec- 
tive distance to the bifurcation point A A* = |A — X* H \, 
and measure from the deterministic temporal series the 
period, T*, of the oscillations when the system is evalu- 
ated at a distance AA* from the deterministic bifurcation 



point. With this information, in Fig. 5 we plot the be- 
havior of t* with the quantity AT* = \T* - T(X H )\, 
where T(Xjj) is the period of the deterministic system at 
precisely the Hopf bifurcation point. The curves can be 
fitted by a power law r* ~ (AT*) Q with the exponents 
a = -0.59, -0.58 and -0.53 for (g), (g) and (g) respec- 
tively, and this seems to indicate that the localization 
behavior with r is characterized by a unique exponent 
with value close to —1/2. We note that for the case of 
(|l|) it is even possible to relate AT* with the system im- 
plicit periodicity T2, thus recovering a similar relation to 
that calculated in p2| . In this way, these results relate 
the resonantlike behavior previously reported in fl2]| us- 
ing the quality factor (3 [|l9| , with an increase (in mean) 
of the localization of the orbits of the system. 
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FIG. 4. h/W versus r for different (A) for (a) Eq. 
(0) (h/W x 10 3 , o r = 0.05 and, from top to bot- 
tom (A) = 2.00828,2.01055,2.01344,2.01713,2.02182 and 
2.02781), (b) Eq. (g) (h/W x 10 6 , a = 5 x 10~ 4 
and, (A) = 0.124,0.1235,0.123,0.1225,0.122 and 0.1215), 
and (c) Eq. (g) (h/W x 10 5 , a = 10" 3 and, 
(A) = 0.496,0.494,0.492,0.49,0.488 and 0.486). Number of 
realizations, iterations and time steps are as in Fig. 3. 

From the behavior of the quantity h/W, it is clear that 
a concentration of orbits around a narrow range of bands 
in the phase space implies a bigger weight of those partic- 
ular frequencies in the power spectrum, and, as a conse- 
quence, a nonmonotonous behavior qualitatively similar 



to that of Fig. 4 should be expected for /3, indicating an 
increase of the coherence in the system response. This 
is indeed the case for our three models (with quality fac- 
tors showing a maximum for values of the correlation 
time close to r*) clearly indicating that the SR-like ef- 
fect induced by colored noise in nonlinear systems with 
limit cycle behavior is quite general. 

Summarizing, we have presented numerical evidence of 
a novel effect of enhanced localization of orbits mediated 
by the correlation time of a multiplicative OU process in 
nonlinear dynamical systems with limit cycle behavior. 
This effect is characterized by a power law with expo- 
nent close to —1/2 for all the models considered in spite 
of their different nonlinearities. This behavior could indi- 
cate the universal character of this phenomenon, but fur- 
ther research is required to clarify this point. This work 
also relates the SR-like behavior previously reported with 
this localization effect. 

We acknowledge financial support from DGESEIC 
(Spain) project PB97-0076. 
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FIG. 5. r* versus AT* for (a) Eq. (Q), (b) Eq. (|), and (c) 
Eq. (M). Simulation parameters are as in Fig. 3. 
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